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ABSTRACT 


Calculations of the reflectivity of water clouds (liquid 
and ice particles) are compared to observations of terrestrial 
clouds in the near infrared. The presentation is divided into 
four parts which may be consulted individually: §3 presents new 
Mie scattering calculations of general interest, || 4-7 compare 
multiple scattering results to cloud observations, $ 8 suggests 
a revision ip the optical constants of ice for \ 3 lJ , and the 
appendix details several methods which substantially reduce the 
work load in multiple .scattering computations. 

Our results indicate that it is possible to use the spectral 
variation of the reflectivity to derive the size of the cloud par- 
ticles and their phase (liquid or solid) as well as the total 
optical depth of the clouds. Typical results show dense cirrus 
clouds to have an optical depth > 10 and to be composed of ice 
particles of mean radius 15-20 M ; the cumulus clouds which were 
analysed showed a more variable, but usually smaller', particle 
size . 

In spectral regions whore the single scattering albedo is 
high it is found that most of the gas absorption ta<es place within 


the clouds rather than above them. 
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Introduction 


Sagan and Pollack (1967), Pollack and Sagan (1968), and 
Hansen and Cheyney (1968, 1969) have pointed out that 
some physical properties of clouds, such as mean particle 
size and optical depth, can be inferred from an analysis 
of the clouds* near infrared reflectivity. Below we 
analyze observations of terrestrial water and ice clouds with 
this expectation in mind . It affords us an opportunity both 
to deduce the properties of these clouds and to check the 
consistency of our approach when applied to clouds we know 
something about. It is especially important to obtain 
tests of the plane-parallel approximation for real atmos- 
pheric cloud systems and the Mie scattering approximation 
for ice crystals . 

After outlining our computational scheme, we summarize 
the theoretical results for single scattering. These illus- 
trate the effects on the scattering of various input oara- 
meters. Subsequently we discuss the infrared observations 
of Blau, et al. and attempt to remove the effects of gaseous 
absorption. In the next section we summarize the absorption 
corrections and relate them to the clouds' temperature and 
to the mode of line formation. Finally, we estimate the 
phase, mean particle size, and optical depth of selected 
clouds from their spectral and angular scattering behavior. 
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.1* 




2 . Computational method 

To obtain the spectral reflectivity of clouds and the 
angular distribution of the scattered light we first com- 
puted the single scattering from a small volume containing a 
representative distribution of particle sizes and then the 
multiple scattering for the entire cloud. The single scattering 
computations were made using the Mie theory in essentially the 
way described by Deirmendjian and Clasen (1962) . However, the 
logarithmic derivatives of certain spherical Bessel functions 
occurring in the series were obtained by proceeding from the 
highest order terms? Kattawar and Plass (1967) have pointed out 
that an upward progression for the recursion relations is basically 
unstable and may cause significant errors for large size parameters. 
The validity of our single scattering solutions was checked by 
making comparisons to published results of Deirmendjian (1964) , 
unpublished company reports by the same author and unpublished 


work of H. Cheyney. 



We employed Deirmendjian ’ s (1964) "Cloud model" for the 
distribution of particle sizes, 

fl 

n(r) cc r exp (-6r/r m ) (1) 

where n(r) is the volume concentration of radius r and r 
is the radius at which the distribution has its maximum. 
Although sample size distributions from actual clouds often 
differ markedly from Deirmendjian* s model, it serves the 
purpose of averaging out most of the large fluctuations which 
occur in the phase function (scattering diagram) for a single 
sphere. We have made preliminary computations to find the 
effect of changing the shape of the size distribution; the 
results for several different distributions indicate that the 
volume extinction (a ext ) the single scattering albedo (o>), 
the asymmetry factor (<( cos 9 )>) and the shape of the phase 
function (outside the region of the glory) depend mainly on 
the mean particle radius for extinction, 

7 = I^ext ( r > *(r)rdr 

CD * ' ' 

l ff ext n ( r > dr 

ana not on the shape of the size distribution. This indi- 
cates that our computations in this paper with a single type 
of size distribution are meaningful, although in cases where 
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a more exact knowledge of the size distribution is available 

it should of course be used. For this paper our computations 

were made with r m = 2 , 4, 8, 16 and 32jj with the integrations 

over particle size extended to the radius r max = 25(j for the 

first four distributions and to r „ - 50u for the case r„ = 

32^. The integration increments were chosen small enough to 

make the phase function smooth (see, e.g., Dave, 1969a, 1969b). 

The values used for r were somewhat arbitrary but that is 

max 

not essential since the purpose of the integration over particle 

sizes was to smooth out single particle effects? we usually 

found *r =- r„. 

m 

Except where otherwise indicated we employed the optical 
constants for water and ice tabulated by Irvine and Pollack 
(1968) ? the computations were made at each wavelength (~ 50) at 
which the authors tabulate the optical constants for the region 
of interest (1.2 £ X £ 3.6jj) . 

The intensity of light multiply scattered by the clouds 
was obtained by using the "double only" computing method, 
described by Hansen (1969a) , which is a variation of van de 
Hulst's (1963) doubling method. If errors ^ 1% are tolerable, 
as is certainly the case for comparison to most observations, 
then the computing time may be greatly speeded up. Some of 
the more useful shortcuts which we have tested are described 


in the appendix. 
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3. Results for single scattering 

The single scattering behavior of an ensemble of aerosols 
is a function of the complex index of refraction, n c ■ n r - in^, 
and the distribution of particle sizes, expressed in units of 
the ratio of the circumference of the particle to the wave- 
length xs x * 2rrr/x* Assuming that the distribution of size 
is given by (1) , we find that the parameter x^, which equals 
2rrr m /x, serves to define the dependence of the scattering 
behavior upon particle size. The results described below are 
of both general interest and of help in understanding the 
infrared properties of terrestrial clouds analyzed in later 
sections. In the calculations for Figs. 1-3 the integrations 
over the size parameter x were performed to an upper limit of 

x - 2x , except when x * 128 for which x was set equal 
max m m max 

to 200. 

a. General single scattering results 

Fig. 1 shows the dependence of the normalized phase func- 
tion p(6) (Deirmendjian, 1964, 1969) upon the particle size 

parameter x in the case of no absorption (n. * 0) with the 
m i 

real refractive index equal to 1.33, the value for liquid water 
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in the visible. The variable 9 is the angle of scatter. 

The individual curves have been vertically displaced 

from one another. The short horizontal line intersecting 

a given phase curve denotes the position at which the phase 

function has a value of unity. For the largest values of 

x ffl , the very precipitous decline in the value of p near 0° 

corresponds to the diffraction peak (Deirraendjian, 1964, 

Dave, 1969c) ; the strong maximum at a scattering angle o£ 

142° is the main rainbow (caused by rays undergoing a single 

internal reflection) with its first supernumery bow (van de 

Hulst, 1957, p. 241) located at 147°? the second rainbow 

(2 internal reflections) is located at 123° while its first 

supernumery bow lies at 114°; and finally the glory corresponds 

to the overall increase and oscillatory behavior near 180°. 

As x decreases all these features become less pronounced and 
m 

broader; in addition the rainbow shifts in location towards 
larger scattering angles and the slope of the phase function 
decreases until near x m * 1/2 it is very similar to the Rayleigh 


phase function. 



Ihe effect of Introducing absorption within the particles 
is investigated in Fig. 2. For x m = 32, the glory and rain- 
bow have been effectively suppressed when n^ is comparable 
to or larger than 0.03. This effect can be understood using 
concepts from geometrical optics. A ray traversing a path 
equal to the particle radius will be diminished in strength 
by exp(-ka) = exp (-2 x n^) where k is the linear absorption 
coefficient. Setting n^ equal to 0.03 and x to 32, we see 
that the ray is diminished by almost a factor of 10 in in- 
tensity. Similarly the diffraction component which is not 
affected by absorption, extends to larger scattering angles 
as n^ increases. 

For x m again of 32, the value of the phase function in 
the backward hemisphere (9 > 90°) begins to increase as n^ 
becomes larger than 0.1. This may be attributed to an en- 
hanced value of externally reflected light, as seen from the 
Fresnel equation. For the smaller particles (x^ = 2) the 
results are basically similar, although less pronounced. 

Finally we study the influence of the real part of the 
index of refraction in Fig. 3. For x m = 32, increasing n p 
causes the rainbow to shift to larger scattering angles, 
until it merges with the glory, resulting in a large increase 
in p (180°). Aside from the diffraction peak, small angle 
scattering becomes more dominant as the refractive index 
decreases, a result which shows up in the asymmetry factor 

< cos 9 > described below. We will find this effect to 



be of significance in understanding certain spectral charac- 
teristics of water clouds. Similarly at the lover values of 
n r refraction tends to dominate over diffraction a t smaller 
angles of incidence and hence the break in the diffraction 
peak occurs at smaller values of 0. The right hand half of the 
figure shows that similar effects occur for particles with 
absorption (n^ * 0.01) and small particles (x m * 2). 
b. Single scattering for water and ice 

One important integral scattering parameter is the single 
scattering albedo, u;, the ratio of the amount of light scattered 
to that which is scattered and absorbed. In Figs. 4 and 5 the 
spectral variation of 1 - w is exhibited for liquid water and ice 
particles, respectively, for five values of the size parameter r m . 

As noted by Irvine and Pollack (1968) , the maxima and minima 
of ice and water are displaced by about O.i^ in wavelength, a 
feature useful in distinguishing the phase of water clouds. Also 
the single scattering albedo generally declines monotonica.1 ly with 
increasing particle size* It should be noted that in all cases 
x m > 3 . For sizes much smaller than this the particles will become 
completely absorbing and the above generalization is no longer true 
(van de Hulst, 1957) . An interesting mild deviation from the 
general variation of £ with particle size occurs near a wavelength 
of 3|a and may be due to the perturbing effect of the smallest 
particles having values of x <1. Near 3p the value of n^ is so 
large that the particles become completely opaque so that allowing 


for the diffraction peak Z is about 1/2. 

The asymmetry factor <( cos 9 )> is defined ar the solid 
angle average of cos 0 weighted by the phase function. It 
describes the degree of forward scattering of the phase 
function. For isotropic scattering it has a value of zero 
while it approaches unity as small angle scattering tends to 
dominate. Figs. 6 and 7 show the spectral behavior of 
< cos 0 y for liquid water and ice particles, respectively. 
For the largest value of ^ there is a very pronounced 
increase in <( cos 0 y slightly shortward of 3 m- and a less 
obvious minimum somewhat longward of 3^. This behavior 
reflects the anamolous despersion changes in n r near the 
very strong absorption feature centered near 3p, with change* 
in also influencing <f cos 0 > nea~ its minima. As shown 
in Fig. 3 and discussed earlier the phase function becomes 
more forward scattering as n r approaches unity. 

From multiple scattering computations reported below 
and elsewhere (Hansen, 1969 a; Hansen and Cheyney, 1968 ) 
we find that for wavelengths less than 2.5u and between 3.4 
and 3.6u the reflectivity of thick clouds depends primarily 
on the single scattering albedo* qualitatively the curves 
of log (l-o)) are very similar to the spectral variation 
of the cloud reflectivities. As the single scattering albedo 
decreases, the less probable it is for a photon to survive 
a number of scattering events. Since the single scattering 



albedo in these wavelength regions varies systematically 
with the characteristic particle size r m (cf. Figs. 4 and 3) 
we can obtain particle size information from the near infra- 
red spectral behavior of clouds, as has been pointed out by 
Sagan and Pollack ( 1967 ), Hansen and Cheyney (1968), and 
Irvine and Pollack ( 1968 ) . 

On the other hand in the spectral region between 2.9 
and the single scattering albedo has a constant value 

of about 1/2 and the spectral behavior of the cloud re- 
flectivity depends primarily on < cos 9 )>. There is an 
inverse qualitative similarity between the wavelength varia- 
tion of < cos 9 > and computed cloud reflectivities. As 
<( cos 9 y increases, less light incident on a cloud' layer is 
reflected back out in the first few scattering events . If 
the single scattering albedo is sufficiently 3e ss than one 
so that these scattering events are the major producers of 
reflected photons, the cloud reflectivity will decrease as 
<( cos 9 y increases . For example the maximum in < cos 0 > 

slightly shortward of 3q results in a minimum value for the 
cloud reflectivity at the same wavelength, as shown below. 
The spectral variation of < cos 0 > is a reflection of the 
spectral behavior of the indices of refraction. 
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4. Reflectivity measurements of terrestrial clouds 

Blau, Esp.inola, and Reifenstein ( 1966 ; see Blau and 
Espinola, 1965, for a more detailed report) have obtained 
airborne infrared measurements of the reflectance properties 
of terrestrial clouds. These consist of both spectral and 
angular observations. Observations were performed from an 
aircraft above the cloud of interest and on a given day- 
measurements were made either in the 1.2 p to 2.5m- region or 
the 2.4p to 3.6p, wavelength domain. In addition on some 
occasions angular scattering information was obtained by 
flying along a hexagonal path and performing spectral measure 
ments of the cloud area situated at the center of the hexa- 
gon. In such measurements the angles of refj.ection and 
incidence remain constant, while the azimuth and angle of 
scatter vary. Below we describe these observations in greate 
detail and outline the transformation we applied to them so 
as to be able to compare them with our multiple scattering 
computations . 

The spectral radiance (specific intensity) 1^ values 
reported by Blau et al. refer to averages of a number of 
spectra obtained close together in time. In addition they 
also give values for the standard deviation of each average 
value. The standard deviation is not necessarily a reflec- 
tion of the error of measurement, because it also includes 
short term variations in the properties of the observed cloud 
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We subjected these measurements to two types of transforma- 
tions. First we divided the observed specific intensity by 

the solar flux outside the atmosphere of the earth, F.: 

A 


R 


A 


^O X/ 7T 


s 

^O 


( 3 ) 


where S, defined by Eq. (3) , is the usual scattering function, 

cos 1 u o the angle of incidence, and 0 = cos _1 fj the angle of 
reflection. If the clouds were a lambert surface, i.e., 
if they scattered isotropically, the normalized reflectivity 
would equal the spherical albedo of the clouds and would 
exhibit no angular variation. 

All the spectra show absorption features due to water 

vapor. In addition some absorption due to carbon dioxide 

at is expected and there are strong C0 2 absorption features 

near 2.7 (i. Because gaseous absorption takes place not only 

above the clouds but through multiple reflection within the 

clouds, we did not exactly correct for this effect. One 

cannot practically correct for the multiple scattering within 

the clouds by assigning an effective single scattering albedo 

to the gaseous absorption component, because the spectral 

resolution was much larger than the spectral domain over 

which the gaseous absorption is constant. We performed our 

water vapor absorption corrections by comparing the value of 

R at two spectral positions, one position expected to have 
A 



very little gaseous absorption and the other a large amount. 
Furthermore the clouds were expected to have approximately 
the same intrinsic value for R^ at the two wavelength posi- 
tions. In the 1.2 - 2.5 p, region we compared values at 
1.28(i with either 1.35p or 1.58p., while for the 2,4 - 3.611 
domain we considered values either at 2.50 or 2,55 with ones 
at either 2.60 or 2.6l(i. The ratio of the values of R-^ were 
equated to ratios of gaseous transmissivities calculated by 
Wyatt, Stull and Plass (1962). The comparison was made at 
pressures and temperatures closest to the cloud top condi- 
tions, as inferred from the cloudtop altitude, and in all 
cases at an effective resolution of lOOcm”^. While this 
resolution is somewhat poorer than the resolution of the 
spectrometer it was found to yield the most consistent 
results . In part the need to employ 100cm resolution 
is a reflection of the breakdown of the theoretical statistical 
model at finer resolutions, as indicated by the appearance 
of high frequency features whose amplitude is too large. 

From the comparison with Wyatt et al.’s tables we derived 
an effective water vapor abundance, W, and used this amount 
to estimate the effective gaseous transmission at other 
wavelengths. At wavelengths where there is more absorption 
by the cloud aerosols, there is less multiple scattering and 
in this sense the transmission correction is an overestimate. 
This circumstance pertains at most of the other wavelengths 
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and the true cloud reflectivity in general should lie between 
the uncorrected and "corrected" values. 

* Absorption effects by carbon dioxide were much more 
localized in the spectrum. For the strong absorption band 
near 2.7 m,, we derived an effective amount of C0 2 and an 
average pressure by allowing for the path down to and up from 
the cloud top and estimating the effective path in the clouds 
from the water vapor absorption amount. The latter can be 
related to the atmospheric temperature where the absorption 
takes place, as discussed below, and to the pathlength with 
the help of the U.S. Standard Atmosphere Supplements (19 66). 

For the much smaller corrections within the 2[i band, we used 
the water vapor absorption amounts as a guide in a less rigor- 
ous fashion. Transmission corrections were then obtained 
from the tables of Stull, Wyatt and Plass ( 1963 ) . 

Fortunately the angular measurements refer to the path- 
lengths of nearly constant angle of incidence and reflection 
and so to first order they require no absorption correction. 
However because some absorption takes place through multiple 
scattering the effective atmosphere transmission may have an 
azimuthal dependence. This effect is very difficult to correct 
for and no attempt was made to do so . 
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5. Discussion of gaseous absorption corrections 

In Table 1 we summarize the water vapor amounts, W eff , 
deduced from Blau et al.’s spectra. The first column gives 
the figure number of the spectra in their final report. 

P is the pressure assumed in deriving W ef . f . Also given are 
the cloud type and the cloudtop altitude. Using the U.S. 
Standard Atmosphere Supplements (19 66), we have estimated 
the temperature at the cloudtop, T Qd , from its altitude. 

To assess the contribution of multiple scattering within 
the clouds to the observed water vapor absorption features, 
we have computed the equivalent amount of water vapor W ^ 
which the sunlight passes through above the clouds, on its 
path down to the cloudtop and up to the airplane. In per- 
forming the calculation we have assumed the atmosphere to be 
saturated. In addition we have corrected for the difference 
in the value between the pressure used in obtaining W ^ 
and the actual cloud top pressure by assuming that pressure 
and gas amount ai*e equally effective in causing absorption. 
For the first three spectra of the 1.2 - 2.5|i region W ref 
is significantly smaller than . Hence most of the absorp 
tion takes place in the clouds. A confirmation of this de- 
duction is obtained by comparing the values of W eff for the 
two sets of spectra from Fig. 55. The value of w eff obtained 
for the larger angle of incidence is smaller than that for the 
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smaller angle of incidence; this finding is opposite to the 
expectation for absorption taking place above the cloud, but 
in accord with predictions for multiple scattering within a 
cloud layer (Chamberlain, 1965) . Such a phenomenon is also 
present for at least some of the gaseous absorption features 
of Venus (Chamberlain and Kuipe.r, 1956), 

On the other hand, W f is larger or comparable to W eff 
for the two spectra pertaining to the 2.4 - 5.6 p region. 

This results in part from the cloudtops being located high 
within the stratosphere where the water vapor abundance is 
only a few percent of the saturation abundance . For these 
spectra we can conckide that the fraction of the cloud signi- 
ficantly contributing to the scattering lies within the 
stratosphere and therefore we are surely viewing ice 
particles . The relatively small depth of penetration for 
these spectra is a result of the highly absorbing nature 
of ice aerosols at 2.5 and 2.6ji, the wavelengths at which 


W 


eff 


was obtained . 


To derive an estimate of the depth of penetration for 


the other clouds, we have computed T scat , the base tempera 
ture required within a saturated atmosphere so that light 


traveling on a straight line down to and up from this level 
at a 60° angle would experience the observed amount of absorp 


tion. Since the actual pathlength within the clouds is more 


tortuous, is probably a slight overestimate of the 
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level of line formation. We see that for the first three 
spectra T scat is substantially larger than T implying 
substantial penetration within the clouds. At the wave- 
lengths used to derive W e ^ for these spectra the cloud 
aerosols are essentially transparent. 

6. Analysis of the cloud spectra 

Our theoretical spectra are functions of three parameters 
the characteristic particle size r m , the optical thickness t 
of the clouds, and the phase of the cloud, i.e., whether the 
aerosols are liquid water or ice , Below we attempt to 
estimate each of these parameters by comparing the theoreti- 
cal and observed spectra. 

In Figs. 8 and 9 Blau et al's observations of a cirrus 
cloud in the 1.2 to 2.5u wavelength region are compared with 
theoretical spectra for an ice cloud. Hie observations 
correspond to Blau et al.’s Fig. 55 for an angle of incidence 
of 71°. Circles and solid bars represent the average values 
and standard deviations of the reflectivity R (eq. 5) after 
correction for gas absorption, while the uncorrected observa- 
tions are indicated by triangles and dotted bars in the cases 
where they differ significantly from the corrected values. 

As explained above, the absorption corrections may be over- 
estimates in regions where the cloud aerosols strongly absorb. 
In this case the true reflectivity will be somewhere between 
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the corrected and uncorrected values. The theoretical curves 
in each portion of the diagrams correspond to various choices 
for the optical depth* the particle size parameter is vailed 
by a factor of 2 between successive sections. Comparing the 
four sections of Figs. 8 and 9, we see that the observations 
permit the determination of a well defined value for r m , l6y.. 

In all of the figures the optical depth refers to "K = 1.2 p.,* 
however, the wavelength dependence of t is small since we 
always consider a distribution of particle sizes with x m >l. 

In determining t, values of R between 1.2 and 1.% are of 
particular use: the theoretical curves differ at these 

wavelengths by a maximum amount. For the cloud observations 
illustrated in Fig. 9 a value t ~ 10 appears to give the 
best fit. However, the optical depth is more difficult to 
estimate than the particle size, and this derived optical 
depth should probably only be regarded as a lower limit. 
Danielson, Moore and van de Hulst (1969) have argued, from 
observations and computations, that condensation nuclei 

limit the cloud reflectivity at wavelengths where ice and 

_ 

water do not effectively absorb [1 - o> <,10 ' ]• 

As mentioned above the absorption features of water and 
ice are displaced somewhat from one another. For example 
peak absorptions and hence minimum values of R occur at 1.45 
and 1.95u for water and at 1.52 and 2.00p for ice. As a 
result the theoretical water spectra do not fit the observed 



20 


spectra in these regions. On the other hand, there is some 
indication that a local maximum in R occurs at 2.2^, a posi- 
tion expected for water clouds, rather than 2.3jj, the place 
for ice clouds. A similar situation was found in some of the 
other spectra. This could be understood in terms of a mixed 
phase model with ice predominating near the top of the cloud, 
but the quality of the present data did not warrant calculations 
for such a model. 

In a similar fashion we analyzed several other observed 
spectra. The results are summarized in Table 2. It is en- 
couraging to see that the same cloud properties were found from 
observations of the same cloud viewed at two different solar 
illumination angles (Fig. 55) . An interesting feature of 

Table 2 is the similar value for r M found for the various ice 

m 

cloud spectra. The deduced average particle radius of 16 M , or 
diameter of 32p, is consistent with values typically obtained 
from direct sampling measurements of ice clouds. 
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TABLE 2 

DEDUCED PROPERTIES OF SOME OBSERVED TERRESTRIAL CLOUDS 


Blau ' s 
Fig. 


Cloud 

Deduced Properties 
Phase r m (u) t 



A = 1,2 to 2. 




55 

71 

cirrus 

ice> 

16 

10 

55 

64 

cirrus 

ice 

16 

10 

53 

44 

cumulus 

water 

8 

4 

51 

52 

thick cirrus 

ice 

16 

10 



over cumulus 






X = 2.4 to 3. 

6u 



18 

48 

cumulus 
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In Fig, 10 we consider the degree to which multiple 
scattering occurs in the 2.i to >.5u spectral region, 
where the cloud aerosols are highly absorbing. The curves 
illustrated are for an optical depth of unity and a large 
enough optical depth (t = 128) to be considered equivalent to 
an infinite value ; we obtained almost identical results for 
all optical depths in excess of 10. The dashed curve repre- 
sents photons scattered only once, the dotted curve photons 
scattered n times (n> 1) with the first n-1 scattering events 
ccc-uring withui the diffraction peak (Hansen, 1969 b), while 
the solid curve represents all the photons reflected from 
the cloud . We see that even if the dotted curve is considered 
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as representing single scattering, single scattering computations are 
inadequate to describe the reflectivity in this highly opaque part of 
the spectrum. 

In Fig. 11 we compare computed cloud spectra for t > 10 with data 
given in Blau et al.*s Figs. 18 and 20. The theoretical curves corre- 
spond to various choices of r m . The data points between 2.5 and 2.6 p 
indicate a particle size of about 16^ in both cases. Near 3fi the ob- 
servations and calculations are not in good agreement; this disagreement 
is discussed in detail in $ 8. 

Finally we note that sample calculations of blackbody thermal emis- 
sion from the clouds showed this contribution to be very small compared 
to reflected sunlight, even at the longest wavelengths of observation. 

7. Analysis of the clouds' angular scattering 

In this section we attempt to assess the information content of 
the clouds' angular scattering properties. As mentioned earlier, ob- 
servations were made at nearly constant angles of incidence and re- 
flection but with varying azimuth and hence scattering angles. In 
Figs. 12 and 13 we contrast the theoretical single scattering behavior 
of water clouds with the complete multiple scattering behavior. The 
computations were made for angles of incidence and reflection of 60° 
and 80°, respectively. In Fig. 12 the reflectivity has been calculated 
for a distribution of water particles having a mean size r m of 2p, while 
Fig. 13 pertains to a value of 16 M for r m . For the larger size particles, 
we defined single scattering to include photons scattered (n-1) times 
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in the diffraction peak (n ;> 1) before suffering a final 
scattering back out of the cloud. Since the 2\i size 
particles do not have a striking diffraction peak, we 
defined single scattering for them in the conventional 
manner . 

We see that the effect of multiple scattering is to wash 
out features such as the rainbow and to greatly diminish the 
angular variation of R. Even for an optical depth of unity 
the actual scattering behavior is markedly different from 
that of single scattering. At an optical depth of 16 the 
scattering behavior is quite close to that for an infinite 
optical depth. On the basis of their high albedo in the 
visible we know most clouds have an optical depth of 16 or 
more. In this event a knowledge of the exact value of the 
optical depth is not too important for being able to predict 
the reflectivity properties . 

In Fig. 14 we compare theoretical calculations with the 
angular dependence exhibited by a cumulus cloud whose top 
was at 4 kft. This doud was part of Hurricane Gladys. For 
each data point the angle of reflection was 80° but the 
angle of incidence varied between 56° and 64°. The 
theoretical computations were performed for spherical water 
particles at the angles appropriate for each observation, 
and the theoretical points were joined together by straight 
lines. As mentioned above, no correction for gas absorption 


was made 
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All the theoretical spectra are normalized to fit the 
data point at the lowest value of (0 - P c ) • A particle size 
of B^l appears to yield the best fit to the measurements, a 
result compatible with typical direct sampling measurements 
of cumulus clouds. Unfortunately, the absorption corrections 
were so large we were unable to meaningfully analyze the 
spectra to confirm this deduction. 

Angular computations in the 2.4 to 3.6p, region are shown 
in Figs. 15 and 16 . The calculations were made for 
spherical ice particles with t he realization that this may 
lead to a very poor approximation to the true angular scatter- 
ing behavior of ice clouds. In Fig, 15 a mean particle 
radius r m of 16g, wasused and the optical depth varied. 

In Fig. 16 the particle size is varied, while the optical 
depth was set equal to 32. In all cases the angles of 
incidence and reflection are 50° and 80° respectively. 

We see that for cases of intermediate aerosol absorption, 
e.g., r m = l6u and ~h = 2.J2\i, the exact solutions most 
clearly preserve such features as the rainbow peak and 
exhibit the strongest dependence upon particle size. When 
the aerosol absorption is relatively small, many scattering 
events occur leading to a smooth, featureless angular 
behavior. When the aerosol absorption Is very large the 
rainbow is not present even for single scattering. 

In Fig. 17 we compare theoretical and observed cloud 
reflectivities for an angle of reflection of 80° and an 
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angle of incidence between 44° and 56°. The observations pertain to cir- 
roform clouds whose tops were at 43 kft. They were part of Hurricane 
Gladys. The theoretical curves have the same meaning as those in Fig. 14, 
except that Fig. 17 is for spherical ice particles. A particle radius of 
16>i i a value found typical for cirrus clouds from the spectral analysis, 
gives no worse a fit to the data than other sizes. The calculations at 
all 3 wavelengths were made with the optical constants given by Irvine 
and Pollack? however, if we accept the modification in the optical constants 
which is indicated by the results in J 8, then the theoretical calculations 
at >. - 3.10(j would be changed and brought into better agreement with the 
observations. At 3.10p the primary effect of the new optical constants 
is to increase the Fresnel reflection in the backward direction; this 
would increase the theoretical reflectivity at the two right-most data 
points in Fig. 17. 

8. Optical constants of ice near \ = 3^ 

In this section we reconsider the ice cloud spectra near \ « 3p where 
the calculations and observations are not in good agreement, and we sug- 
gest one possible explanation for the discrepancy. 

A minimum in the observed spectra (Fig. 11) occurs at about 2.9 jj in 
wavelength and a maximum around 3. Ip. These features are particularly 
prominent in the right hand graph. This may at first appear surprising 
since ice has its maximum absorption at 3.1y (Irvine and Pollack, 1968) . 

Blau and Espinola (1969) first pointed out the minima in the spectra and 
correctly attributed it to the anomalous dispersion of n r , which is dis- 
cussed above. A local maxima in the reflectivity is expected near 3.1(j 
for similar reasons. Our computed spectra qualitatively show these effects, 
but particularly for the right hand graph (Fig. 11) they fail to quantita- 
tively match the observations. Blau and Espinola encountered a similar 
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difficulty in exolainina the minimum and suggested that Irvine and 


Pollack’s value of n r be revised near this wavelength position. 

We have also considered the possibility of revising 
Irvine and Pollack’s optical constants for ice in the vicinity 
of the strong 3*lp- feature. We have attempted to do this in 
a consistent fashion by 'relating changes in the real part of 
the index of refraction, n r , to changes in the imaginary 
part n^ . Spitzer and Kj.einman (1961) have given relation- 
ships between n r and n^ under the assumption that individual 
absorption bands may be considered as classical oscillators 


and have obtained a very good fit to reflection measurements 
of strong infrared bands of quartz. Assuming that the 3.1 jjl 
band is a dominant feature at nearby wavelengths and neglect- 
ing small differences between the value of the frequency at 
the desired positions and the central frequency of the band, 
we have simplified Spitzer and Kleinman’s formulas to the 
following general form: 


pr 2 

n r n i = ( A A) 2 +r^ 


W 




A + 


2pr A Tv 

(A>0 2 +r 2 


(5) 


where p is the maximum value of the product of n^jrn, F one 
half the value of the width between half maximum points of 
n r n ± , AA the difference in wavelength between the value of 
interest and the position of the maximum, and A is a constant. 


These formulae give a good fit to the water and ice data 
of Irvine and Pollack near the absorption feature. 
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For the ice constants given by these authors p is about 0.42, 
and T about 0.15u. The constant A is readily found by apply- 
ing the second formula to wavelengths significantly short- 
ward of jl; we find A is 1.66. 

In modifying the indices of refraction given by Irvine 
and Pollack we assumed that F w$s unchaiged and varied the 
only remaining free parameter p by various scale factors c. 

We have modified the published data between 2,8j± and J,6\i 
because this is the region dominated by the 3.1 p, band. In 
addition the published indices shortward of 2.8p were based 
on a different set of measurements than those used at 2.8p 
and longward. These first set were checked against other 
measurements and found to be in good agreement. The results 
of two trial modifications are given in Fig. 18 along with 
the original values . In addition the effects of the modified 
indices on the single scattering albedo au and the angular 
asymmetry parameter <( cos © )> are shown. At short wavelengths 
the value c = 2 led to values of n p less than 1 and so this 
portion o' the curve has not been drawn. 

Since the optical parameters have not been changed 
shortward of 2.8(jl, the deduced characteristic particle sizes 
and the lower limit on the optical depth will still hold . 

We were able to deduce these from the data points between 
2.47 u and 2.65u. fig. 19 shows theoretical spectra for Blau's 
Figs. 18 and 20 based on values of c = 1 (the old optical 
constants), 1.6, and 2. We see that a value of c of about 
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1.6 leads to a good fit of the spectra for both Figs. 18 and 20. 

Some justification for our method of varying the ice constants 
is provided by the following: For wavelengths < 2.8p the optical 
constants of Irvine and Pollack were based on several sources of 
data which agreed with one another. However, for \ & 2.8p Irvine 
and Pollack had only one source for the ice absorption coefficients 
and this source gave water absorption coefficients for \ ~ 34 which 
were ~ 60% less than the values obtained by several other experi- 
menters. Thus the proposed revision of the indices of refraction 
seems quite reasonable. 

9 . Summary 

This study indicates that characteristic broad band absorption 
features in light scattered by clouds can be used to help identify 
the scattering material and to determine the particle size and cloud 
optical depth. The results are consistent with the assumption that 
it is adequate to employ the spherical particle approximation in 
calculating the spectral reflectivity. As far as the angular be- 
havior is concerned, not enough accurate observations are available 
for an adequate analysis of the theory. More extensive measurements 
on atmospheric clouds are desirable, as well as carefully controlled 
laboratory observations along the lines of the recent measurements 
by Zander (1968) and Plummer (1969) . 
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APPENDIX 


The following shortcuts in the numerical work have been 
tested and found useful for the multiple scattering computa- 
tions : 

(a) For large particles (x = 2ira/\ > 25) there is a 
sharp diffraction peak in the phase function which nec- 
essitates a large number of terms in the cosine expansion 
of the scattering function and a large number of points in 
the integrations over p. = cos 9. However, in the case of 
conservative scattering we have shown elsewhere (Hansen, 1969b* 
see also Potter, 1969 ) that the photons scattered into the 
forward spike may be approximated as being unscattered by 
truncating the forward peak from the phase function and 
reducing the interaction optical thickness t such that 

T* = (1 - F) T (Al) 

where t* is the optical thickness to be used with the trun- 
cated phase function, F is the fraction of photons scattered 
into the forward peak, 

F = 5 (P - P*) |f* ( A2 ) 

4tt 

and p and p* are the untruncated and truncated phase func- 
tions, respectively. For conservative scattering this approximation 


30 



31 

introduces large errors in the reflected intensities only 
for small total scattering angles (0 ~ 0°); it introduces 
errors of a few percent if the total scattering angle corre- 
sponds to a sharp feature in the phase function (such as the 
glory) and it also introduces errors of a few p^r cent if the 
incident or emergent angle Is near grazing (0, 0 Q , ~ 90°) 
elsewhere the error is $ 1$. 

In the case of nonconservative scattering the single 
scattering albedo, a), must be scaled such that 


1 



because the assumption that photons in the forward peak 
should be treated as unscattered implies that the absorption 
cross section is unchanged and the scattering cross section 
reduced by the factor 1 - F. Several tests of this approxima 
tion were made and it was found to become increasingly 
accurate as a) decreased . Even in the wavelength region ~ ^ 
where often ~ 80$ of the single scattered photons were cut 
off with th a truncating of the diffraction peak, the error 
introduced was < 1$ except for scattering angles ~ 0°, 

(b) In the expansion of the scattering function in 
cosines, we writes 

S(t; h, 0; u Q , 0 Q ) = E S m (r; *i, n 0 ) cos n (0 - 0 Q ) (A*) 

m=0 
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where u Q and p, denote the cosine of the angles of incidence 
and reflection, respectively and 0 Q and 0 are the correspond- 
ing azimuth angles. The number of terms, M, needed to obtain 
an accuracy within for all p. and p. with a strongly 
anisotropic phase function is typically about one hundred. 
However, the range of p, and p, 0 for which the numerical value 
of S m (u, p 0 ) is not negligible decreases as m increases 
until for m ~ M only S(tj p. ~ 0 \x Q ~ 0) is significant. 

If S (t; p., 0 t p. Q , 0 Q ) is to be calculated at N values of 
p. on the interval (0,1) and N values of p, Q then S°(ti p, p q ) 

p 

must be computed at N points and in the integrations over p 

which include S°(tj p, P Q ) as a factor, N points are employed; 

however the number of points required decreases steadily as 

M 

m increases until for S (tj p., p Q ) calculations are only 
needed for one point and only one point is needed in the 
integrations. For a given accuracy specification it is 
easy by numerical testing to find the number of points at 
which S m (T; p., p. Q ) must be calculated and included in the 
integrations. A factor of 2 - 3 in computer time may be 
saved with an introduction of errors ~ 1 - 2 

(c) The strongest azimuthal dependence arises from 
single scattered photons but an analytic expression exists 
for the intensity due to these photons (Hansen, 1969a) and 
hence considerable computing time may be saved by writing 



(A5) 
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S(t;^i ,0hi o , 0 q ) = s sg (t^ J 0;n o ,0 o ) 

+ Z [S m (T,u,M- ) - S® (t)n,0 )] cos m(0 - 0 ) 
m =0 o ss o o 

where the subscript ss labels the contributions of single 

scattering to the S function. Typically the number of terms 

needed in (A5) is ~ 50^ of the number required with (A4) 

to achieve the same accuracy, 

(d) Several additional ways to save computer time, 

which we found by numerical experimentation, can be shown 

to have a firm theoretical basis from work of van de Hulst 

(forthcoming book) . van de Hulst shows that each term 

in the cos m (0 - 0 O ) expansion may be thought of as having 

an effective albedo for single scattering and this albedo 

decreases steadily as m increases. Some consequences are: 

the doubling process for terms with m >1 may initiated at 

_ic «.ps 

an optical thickness ^ 2 rather than 2 for m > 1 

the asymptotic value of the scattering function is obtained 
already at t ~ 8; for m > M/4 the sum of the infinite series 
occurring in the doubling equations may be replaced by the 
value of the first term. These simplifications may easily 
reduce the computing time by a factor ~ 3 . 

If all four of the above methods for reducing the com- 
puter time are employed the total time saving is not the 
product of the factors which each gives alone because there 
is considerable overlap. In the computations for this paper 
(performed on an IBM 260/95) we always employed (d) and for 
tae particle size distributions with r m = l6*i and 32 4 we 





used (a). All of the above methods may be worthwhile for 
slower computers and especially for problems such as line 
formation and the multiple scattering of polarized light. 

Note added in revision: Dave (1970) has independently shown 
that the number of terms required in the Fourier expansion 
depends strongly on p and ja Q [shortcut (b) above] and he has 
presented graphical illustrations of this. Dave also makes use 
of the fact that the effective albedo decreases toward higher 
terms in the Fourier expansion [shortcut (d) ] . 
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Fig. 


Fig. 


Fig. 


1 Single scattering phase functions for a size dis- 
tribution of transparent spherical particles with 
a real refractive index typical of water and ice 
in the near infrared $ the curves show the effect 
of changing the characteristic particle size. 

In Figs. 1-3 the vertical scales apply to the 
uppermost curve on the left side and the scales 
for the other curves may be obtained by multipli- 
cation by a power of 10 such that the horizontal 
bar on each curve occurs at p(0) = 1. 

2 Single scattering phase functions for a size 
distribution of spherical particles showing the 
effect of absorption within the particles for 
large particles (x = 32) and particles of 
moderate size (x m =2) 

3 Single scattering phase functions for a size dis- 
tribution of spherical particles showing the effect 
of changing the real part of the refractive index 
for large particles with no absorption (left), 

for particles of moderate size with no absorption 
(upper right), and for large particles with moderate 
absorption (lower right). 
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Fig. 

Fig. 

Fig. 

Fig. 

Fig. 


Fig. 

Fig. 


4 Single scattering albedo for a cloud of spherical 
water particles for five different particle size 
distributions . 

5 Same as Fig. 4 for spherical ice particles. 

6 Asymmetry factor of the phase function for a cloud 
of spherical water particles for five different 
particle size distributions. 

T Same as Fig. 6 for spherical ice particles. 

8 Theoretical cloud reflectivities for 9=0, 

9 o = 71° and 0 - 0 Q ^ 180° for five cloud optical 
thicknesses and two size distributions of 
spherical ice particles. The circles and solid 
bars represent observations by Blau, et al. of 
cirrus clouds at 38, 000 feet after correction for 
gas absorption? the uncorrected observations are 
indicated by triangles and dotted bars ir. the cases 
where they differ significantly from the "corrected’* 
values . 

9 Same as Fig. 8 for two additional particle size 
distributions . 

10 Theoretical cloud reflectivities for 9=0°, 9 q = 48 
and 0 - 0 Q = l80°, for two cloud optical thicknesses 
for the size distribution of spherical ice parti- 
cles having the characteristic particle size r m = 
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Fig. 


Fig. 

Fig. 

Fig. 


Fig. 

Fig. 


11 Theoretical cloud reflectivities for 7 > 10 for five 
size distributions of spherical ice particles. The 
circles represent observations by Blau, et al . of 
"cumulus" clouds at 50,000 feet after correction for 
gas absorption; the uncorrected observations are in- 
dicated by triangles. 

12 Theoretical cloud reflectivities as a function of 

azimuth angle for 0 = 80° and 0 = 60° for a size dis- 

o 

tribution of spherical water particles with r m = 2 . 

13 Same as Fig. 12 for a size distribution with r^ = 16. 

14 Theoretical and observed cloud reflectivities for 
0 - 80° and 56 s 0 Q s 64 normalized to unity at 

the smallest observed value of 0 - 0 O * The theoreti- 
cal calculations are for water particles at values 
of 0 , 0 Q , and 0 - 0 Q , correct for each observed point 
and connected by straight lines. The computations 
were made for t = 32, but are approximately valid for 
t a 10. The observations by Blau, et al . were made on 
cumulus clouds at 4,000 feet above Hurricane Gladys. 

15 Theoretical cloud reflectivities as a function of 
azimuth angle for 0 = 80° and 0 O = 50° for a size 
distribution of spherical ice particles having r m * 16 M . 

16 Same as Fig. 15 but for four particle size distributions 
and for only one cloud optical thickness (t ■ 32) . 
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Fig. 17 Theoretical and observed cloud reflectivities for 
0 = 80° and 44 r <: 0 QI s 56° normalized to unity at 
the smallest observed value of 0 - 0 Q . 

Fig. 18 The optical constants for ice are shown in the upper 
part of the figure with the solid curve representing 
the data of Irvine and Pollack (1968) ? the other two 
curves were obtained by multiplying Irvines and Pollack's 
n r n^ by the factor c in the interval 2.8 ^ x ^ 3.6. The 
lower part of the figure shows the single scattering 
parameters < cos 8 > and uu for the three sets of optical 
constants. 

Fig. 19 Same as Fig. 11 with the theoretical curves for the 
size distribution having r m = 16 jj; the three curves 


correspond to the three choices of c shown in Fig. 18. 
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